Spin-resolved Andreev levels in hybrid superconductor-semiconductor nanostructures 
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The combination of superconductors and low- 
dimensional conductors embodies a rich, yet 
largely unexplored physics. In this hybrid sys- 
tem, macroscopic properties enforced by super- 
conductivity can be controlled through electri- 
cally tunable microscopic degrees of freedom, in- 
herent to a relatively small number of confined 
electrons [l|48|. Here we consider the prototypical 
case of a quantum dot (QD) coupled strongly to a 
superconductor (5) and weakly to a normal-metal 
{N) tunnel probe. We investigate the magnetic 
properties of the lowest-energy, sub-gap states, 
which are governed by a competition between 
superconducting pairing and Coulomb repulsion 
[2, 3, 9 13]. In a magnetic field, only when the 
ground state is a spin singlet, can the Zeeman 
splitting of the (excited) doublet be revealed by 
tunnel spectroscopy. The splitting is strongly in- 
fluenced by a level-repulsion effect with the con- 
tinuum of quasi-particle states; and it can induce 
a quantum phase transition (QPT) to a spin- 
polarized state. Our experimental results, sup- 
ported by theory, hold relevance for current re- 
search on quantum-information devices (l4L llSl ] 
and Majorana fermions in hybrid nanostructures 
[16 20]. 

QDs coupled to S electrodes can be seen as artificial re- 
alizations of quantum impurity systems, a core problem 
in the physics of strongly correlated electron systems. 
In the QD-S system, superconducting pairing correla- 
tions favour the mixing of evenly-occupied charge states 
of the QD. Under the simplest assumption of a QD with 
a single spin-degenerate orbital level, four eigenstates are 
possible: a spin-degenerate magnetic doublet (| t), | 4-)), 
and two BCS-like singlets (V) = -v*\ ti) + u|0), and 
\+) = u\ ti) +v*\Q)) O, [ill- For odd occupation, the 
magnetic doublet, favoured by the on-site Coulomb in- 
teraction, and the |— ) state, favoured by superconduct- 
ing pairing, compete for being the ground state (GS), 
whereas |-f) always lies at higher energies. This com- 
petition is governed by the interplay of different energy 
scales: the charging energy, U , the superconducting gap, 
A, the QD-S* tunnel coupling, Fg, responsible for quan- 
tum fluctuations in the QD charge, and the energy eo of 



the QD level relative to the Fermi energy of the S elec- 
trode (see Fig. la) 0, H, [§-[l3| • The qualitative phase 
diagram in the lower part of Fig. la shows that for suf- 
ficiently large Ts/U the QD stays always in the singlet 
GS. 

Here we investigate the magnetic properties of the 
lowest-energy states in a S'-QD-iV geometry, where the 
N contact acts as a weakly coupled tunnel probe. In this 
geometry, a measurement of the difi'erential conductance, 
dl/dV, as a function of the source-drain voltage, V, en- 
ables a direct spectroscopy of the density of states (DOS) 
in the QD-S system. In such a measurement, the fiow of 
an electrical current involves single-electron fluctuations 
in the QD charge. At low energy, these charge fluctua- 
tions correspond to elementary excitations between the 
GS and the first excited state (ES), namely between (| t)i 
4-)) and |— ), or vice-versa. As a result, sub-gap reso- 
nances in dl/dV appear at eV = ijCI, where C is the 
excitation energy (see Fig. la) [2l| - |24| . These sub-gap 
resonances, often referred to as Andreev levels, are re- 
sponsible for the Josephson effect in S'-QD-5 devices. 

We used devices based on single InAs/InP core/shell 
nanowires (NWs), where vanadium (gold) was used for 
the S (N) contact ^25]. A device schematic and a rep- 
resentative image are shown in Figs. Ic and Id, re- 
spectively. The fabricated vanadium electrodes showed 
A = 0.55 meV and an in-plane critical magnetic field 
w 2 T (x II NW axis). The QD is naturally formed in 
the NW section between the S and N contacts. We find 
typical U values of a few meV (i.e., U /A « 3 — 10). The 
QD properties are controlled by means of two bottom 
electrodes crossing the NW, labeled as plunger gate and 
^-barrier gate, and a back gate provided by the conduct- 
ing Si substrate. To achieve the asymmetry condition 
Ts ^ Fjv, the S'-barrier gate was positively biased at 
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2 V. We used the plunger gate voltage Vpg to vary 
the charge on the QD, and the back-gate voltage Vbg to 
finely tune the tunnel coupling. 

Figure 2a shows a series of dI/dV{Vpg,V) measure- 
ments for three different Fg. The top row refers to the 
weakest F5. In this case, the spanned Vpg range corre- 
sponds to a horizontal path in the phase diagram that 
goes through the doublet GS region (see right diagram). 
Let us first consider the leftmost plot taken at magnetic 
field B = 0. On the left and right sides of this plot, 
the QD lies deep inside the singlet GS regime. Here the 
doublet ES approaches the superconducting gap edge, 
yielding an Andreev-level energy w A. By moving to- 
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FIG. 1. Andreev levels in a hybrid S-QD-A'^ system 
and device description, (a) (Upper panel) Schematics of a 
S-QD-N device with asymmetric tunnel couplings to the nor- 
mal metal (Au) and superconductor (V) leads, Fjv and Fg, 
respectively. A is the superconducting gap, U is the charg- 
ing energy, is the chemical potential of the i lead, and eo 
is the QD energy level relative to (in the Fs — )■ limit, 
the QD has 0, 1 or 2 electrons for eo > 0, ~U < eo < 0, 
eo < —U, respectively). The sub-gap peaks located at ±^ 
represent the Andreev levels. In tunnel spectroscopy mea- 
surements the alignment of ^jv to an Andreev level yields a 
peak in the differential conductance. (Lower panel) Qualita- 
tive phase diagram [9|, lid . depicting the stability of the 
magnetic doublet {\D) = \ t), | i)) versus that of the BCS 
singlet (| — )). (b) Low-energy excitations of the QD-S* system 
and their expected evolution in a magnetic field, B. Doublet 
GS case (left): | t) is stabilized by B and Andreev levels re- 
lated to the transition | t) \ — ) are observed. Singlet GS 
case (right): at finite B, the excited spin-split states | t) and 
I D give rise to distinct Andreev levels with energy and 
(^l, respectively. Ez = IqIubB is the Zeeman energy, where 
\g\ is the g-factor and /is is the Bohr magneton, (c) Device 
schematic. The N and S leads were made of Ti(2.5 nm)/Au 
(50 nm) and Ti(2.5 nm)/V (45 nm)/Al (5 nm), respectively. 
The QD system is tuned by means of three gates: a plunger 
gate (pg), a barrier gate (sg) close to the S contact, and a 
back gate (bg). B is applied in the {x,y) device plane {x be- 
ing parallel to the NW) . (d) Scanning electron micrograph of 
a S-QD-N device. 



wards the central region, the two sub-gap resonances ap- 
proach each other and cross at the singlet-doublet phase 
boundaries, where C = 0. In the doublet GS regime be- 
tween the two crossings, the sub-gap resonances form a 
loop structure with ^ maximal at the electron-hole sym- 
metry point. Increasing corresponds to an upward 
shift in the phase diagram. The middle row in Fig. 2a 
refers to the case where Tg is just large enough to sta- 
bilize the singlet GS over the full Vpg range (see right 



diagram). At i? = 0, the Andreev levels approach the 
Fermi level without crossing it. A further increase in Ts 
leads to a robust stabilization of the singlet GS (bottom 
row). At zero- field, the sub-gap resonances remain dis- 
tant from each other coming to a minimal separation at 
the electron-hole symmetry point (eo = — C//2). Similarly 
to the case of superconducting single-electron transistors 
[26j . the QD occupation increases with Vpg in units of 
two without going through an intermediate odd state. 

We now turn to the effect of B on the Andreev levels 
(middle and right columns in Fig. 2a). Starting from the 
weak coupling case (top row), a field- induced splitting of 
the sub-gap resonances appears, yet only in correspon- 
dence of a singlet GS. This is due to the fact that these 
resonances involve excitations between states of different 
parity. For a singlet GS, the spin degeneracy of the dou- 
blet ES is lifted by the Zeeman effect resulting in two 
distinct excitations (see Fig. lb). By contrast, for a dou- 
blet GS, no sub-gap resonance stems from the | t) — ^ I i) 
excitation, because these two states have the same (odd) 
number of electrons. The energy of the only visible An- 
dreev level, associated with the | t) ^ \~) transition, 
increases with B. The formation of a loop structure in 
the rightmost panel of the middle row shows that a QPT 
from a singlet to a spin-polarized GS can be induced by 
B when the starting ( is sufficiently small. In the bot- 
tom row, Zeeman-split Andreev levels can be seen all 
over the spanned Vpg range. At Bj; = 0.4 T, the inner 
levels overlap at the Fermi level, indicating a degeneracy 
between the | f) and |— ) states. In Fig. 2b we show 
theoretical dl/dV plots of a S-QD-N Anderson model 
calculated by means of self-consistent Hartree-Fock the- 
ory [131 (see Suppl. Material). The full phenomenology 
explained above is recovered, supporting our interpreta- 
tion in terms of spin-resolved Andreev levels and a B 
field-induced QPT. 

Interestingly, the splitting of Andreev levels appears 
to be gate dependent. It tends to vanish when the sys- 
tem is pushed deep into the singlet GS, and it is maximal 
near the phase boundaries. To further investigate this be- 
haviour, we have measured dI/dV{B, V) for fixed values 
of Vpg. These measurements were carried out on a sec- 
ond similar device (see Suppl. Material). The mid- panel 
of Fig. 3b displays the B^ dependence of the sub-gap 
resonances measured at position 1 in Fig. 3a. Initially, 
the energy of the Andreev levels increases, as expected 
for a doublet GS (see left panel). From a linear fit of 
the low-field regime, i.e. Ci^x) — C(0) + Ez/'i, where 
Ez — IgxlfJ-sBx is the Zeeman energy and /is is the 
Bohr magneton, we obtain a g-factor \gx\ « 5.6. For 
Bx > 0.7 T, the field-induced closing of the gap bends 
the Andreev levels down to zero-energy. Finally, above 
the critical field, a split Kondo resonance is observed, 
from which \gx\ ~ 5.5 is estimated, consistent with the 
value extracted from the Andreev level splitting. The 
mid-panel of Fig. 3c displays a similar measurement 
taken at position 2 in Fig. 3a, where the GS is a singlet. 
The splitting of the Andreev levels is clearly asymmetric. 
The lower energy level decreases to zero according to a 
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FIG. 2. Andreev levels in different coupling regimes and their magnetic-field dependence, (a) Experimental plots 
of dl/dV vs. {Vpg,V) for different QD-S couplings, Vs (increasing from top to bottom), and different B values (increasing 
from left to right). Top-left panel: along the Vpg range, the system GS changes from singlet (| — )) to doublet {\D)) and back to 
singlet, following the green trajectory in the qualitative diagram on the right side of the same row. We find that increasing Vhg 
results in larger Ts, thereby leading to an upward shift in the phase diagram. Eventually, the green trajectory is pushed into the 
singlet region (mid and bottom diagrams). Experimentally, this results in the disappearance of the doublet GS loop structure, 
as shown in the mid-left and bottom-left panels. The middle and right columns show the _B-evolution of the Andreev levels in 
the three coupling regimes. For relatively weak coupling (top row), the Andreev levels for a singlet GS split due to the Zeeman 
effect, whereas those for a doublet GS simply move apart. At intermediate coupling (middle row), B induces a quantum-phase 
transition from a singlet to a spin-polarized GS, as denoted by the appearance of a loop structure (right panel). At the largest 
coupling (bottom row), the Zeeman splitting of the Andreev levels is clearly visible all over the Vpg range. The splitting is 
gate-dependent with a maximum in the central region, (b) Theoretical dl/dV plots calculated by means of a self-consistent 
Hartree-Fock theory. In all calculations we used U = 2.5A and Fiv/rs — 1/3. From top to bottom, Ts was set to 0.2A, 0.7A 
and 0.9A. In spite of the relative simplicity of the Anderson model, the full experimental phenomenology is recovered. 



linear dependence (f{Bx) = ({0} — Ez/2, with \gx\ ~ 
6.1, which is close to the value measured from the split 
Kondo resonance in the normal state. The higher energy 
level, however, exhibits a much weaker field dependence. 
Both the non- linear field dependence for > 0.7 T in 
Fig. 3b and the asymmetric splitting in Fig. 3c can be 
explained in terms of a level-repulsion effect between the 
Andreev levels and the continuum of quasiparticle states. 
This interpretation is confirmed by numerical calcula- 
tions shown in the right panels of Figs. 3b and 3c, which 
are in good agreement with the respective experimental 
data. In the mid-panel of Fig. 3c, the inner sub-gap 
resonances cross around 1.5 T, denoting a field- induced 
QPT. Above this field, however, they remain pinned as 
a zero-bias peak (ZBP) up to ss 2 T. This peculiar 
behavior, reproduced by the numerical results, can be at- 
tributed to the level-repulsion effect discussed above, in 
combination with the rapid shrinking of A with B^. 

In order to observe a clear _B-induced QPT from a sin- 
glet to the spin-polarized GS, we reduced C(0) by tuning 
Vpg closer to the singlet-doublet crossing in Fig. 3a. The 
corresponding data are shown in Fig. 4a. Contrary to 
the case of Fig. 3c, the Andreev level splitting is rather 
symmetric, owing to the reduced importance of the level 



repulsion effect at energies far from A. The inner sub- 
gap resonances split again after the QPT, which occurs 
now at Bx « 0.5 T. As expected, the outer sub-gap res- 
onances get simultaneously suppressed (left panel of Fig. 
3c). The suppression is not complete though, suggesting 
a partial population of the |— ) ES, possibly favored by 
thermal activation. 

We note that the ZBP at the QPT appears to extend 
on a sizable field range AB^ « 150 mT. This range is con- 
sistent with the r TV-dominated lifetime broadening of the 
Andreev levels, i.e. Ig^-l^sASj; « peak width « 50^eV. 
In Fig. 4b we show how the ZBP depends on the in-plane 
B angle, 9, relative to the NW axis. As 9 varies from 
to 7r/2, the ZBP splits into two peaks with decreasing 
height. This angle dependence is an effect of the g- factor 
anisotropy. For 9 ~ 7r/2, we find a y-factor \gy\ w 3, i.e. 
a factor of 2 smaller than for 6* = (see Suppl. Material) . 
As a result, the QPT only occurs at a higher field (see 
Suppl. Material, Bgpj, sa 1 T), and the split peaks cor- 
respond to Ct transitions on the singlct-GS side. Figure 
4b shows also a pair of small outer peaks associated with 
the Ci transitions. Their oscillatory position is as well 
due to g-factor anisotropy. 

It is interesting to note that the B dependences dis- 
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FIG. 3. Magnetic-field evolution of the Andreev levels at fixed gate voltage and the level repulsion effect, (a) 

dl /dV{Vpg,V) measurement at i? = corresponding to a singlet-doublet-singiet sweep. (b) Left panel; Qualitative _B-evolution 
of tiie low-energy states of a QD-S system as expected for a doublet GS. Middle panel: Corresponding experimental data 
measured at position 1 in (a). ( increases linearly with B until it approaches the edge of the superconducting gap. The 
levels then move towards zero following the B suppression of A. Right panel: Corresponding theory plot of dl /dV{B, eV/ A) 
calculated using experimentally measured parameters, i.e., U/A = 8, eo/A = —4, and Fs/A = 1.5 (Fiv/Fs was set to 1/3). 
(c) same as (b), but for singlet GS. The experimental plot in the middle panel was taken at position 2 in (a). It shows an 
asymmetric splitting of the C^f and peaks. The weak B dependence of is due to the level repulsion between | 4-) and the 
continuum of quasiparticle states above A. The corresponding numerical results were taken with the same parameters as in 
(b), except eo/A = —6. 



cussed above mimic some of the signatures expected 
for Major ana fermions in hybrid devices [H, IT7l [28l - l30| . 
First, a ZBP is observed at finite B along the NW ex- 
tending over a sizable field range. The peak splitting 
beyond the QPT resembles the energy splitting of Majo- 
rana states expected to occur for an increased coupling 
between neighbouring MFs resulting from the extended 
coherence length at higher fields |28l - l3G |. Finally, the 
ZBP is suppressed for 9 = 7r/2, where B is pres umably 
aligned to the Rashba spin- orbit field, S50 [11 El- A 
closer inspection shows, however, that some of the QPT 
signatures differ from those expected from MFs. While 
Majorana ZBPs should occur only for fields perpendicu- 
lar to Bso, QPT- related ZBPs can be observed for any 
field direction, provided the field magnitude is adjusted 
to counterbalance g-factor anisotropies. Secondly, Ma- 
jorana ZBPs are predicted to extend over larger field 
ranges, whereas the field extent of a QPT is limited by 
the ratio between the Andreev-level linewidth (which can 
be as low as the thermal energy in the small-Fjv limit) 
and the g-factor. Figure 3b, however, shows that this lat- 
ter limit can be overcome when the QPT occurs in con- 
comitance with a closing superconducting gap. A similar 
stretching of the ZBP may also arise from nearby sub- 
gap states producing a level-repulsion effect on the inner 



Andreev levels. 



Methods 

Device fabrication. The S'-QD-A^ devices used in 
this study were based on individual InAs/InP core/shell 
nanowires grown by thermal evaporation (3l| (diameter 
w 30 nm, shell thickness « 2 nm). The NWs were de- 
posited onto Si/Si02 substrates on which arrays of thin 
metallic striplines [Ti(2.5 nm)/Au(15 nm), width « 50 
nm] covered by a 8 nm-thick atomic layer deposition 
(ALD) Hf02 film had been previously processed. During 
the measurements, the degenerately-doped Si substrate 
was used as a global backgate, whereas the striplines were 
used as local gates. Normal metal [Ti(2.5 nm)/Au(50 
nm)] and superconductor [Ti(2.5 nm)/V(45 nm)/Al(5 
nm)] leads were incorporated to the devices by means 
of standard e-beam lithography techniques (lateral sepa- 
ration w 200 nm). 
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FIG. 4. Magnetic-field induced QPT and angle 
anisotropy. (a) Left panel: dl /dV{B,V) taken at the po- 
sition of the vertical line in the inset (same device as in Fig. 
3). Right panel: line traces at equally spaced B values as ex- 
tracted from the data in the left panel (each shifted by 0.005 
x2e^/h). The QPT induced by the field is observed as a ZBP 
extending over a B range of about 150 mT. This apparently 
large extension is a consequence of the finite width of the An- 
dreev levels, (b) dI/dV{V) traces taken with \B\ = 0.6 T, 
at different angles. This field magnitude corresponds to the 
QPT field when B is ahgned to the NW axis at 6* = Owing 
to the p-factor anisotropy, the ZBP associated with the QPT 
is split and suppressed when B is rotated away from the NW 
axis. The peak splitting has a maximum a.t 6 = -n /2 
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' Model and Hartree-Fock theory 

D ■ 



P^: . . . 

^ , As described in the main text, our nanowires are well described by a single- level quantum dot with competing 
I • superconducting correlations and on-site Coulomb interactions. A minimal model for describing this experimental 
configuration is giyen by a single Icyel Anderson model coupled to a superconducting reservoir with Hamiltonian 
' 1^- Hqd-s — Hs + -I- Hqe). Here, Hqo models the uncoupled quantum dot, and is given by 

■ Hqd = ^ eadld^ + Un^ni, (1) 



c/2 

CD 



o 



where dj. creates an electron with spin a on the dot level located at e„ . In the presence of an external magnetic field 
B, spin-degeneracy is broken and the levels are given by = eg — ^Ez and = eo + ^Ez with Ez = gfJ-sB, being the 
Zeeman energy. The second term in Eq. ([Ij describes the local Coulomb interaction for two electrons with opposite 
spin within the dot {ua = dlca), where U is the charging energy. Hs describes the uncoupled superconducting lead, 
^ • modelled by a BCS Hamiltonian 

fescr fes 



where = ^fes ~ Ms is referred with respect to the chemical potential at the superconducting reservoir and A is 
the superconducting order parameter. describes the coupling between the QD level and the superconductor and 
has the form 



ks<y 



psj _ This coupling to the superconducting lead is characterized by the parameter Ts — '"'X^fes l^fcsP'^('^ ~ ^fes)- 

■ described in the main text, the competition between the three energy scales U, A and Ts governs the ground state of 

' the model which can be either a singlet or a doublet. Finally, the experimental setup includes a normal reservoir. We 

I • model this by adding two extra terms to the Hamiltonian such that the total model reads H = Hn + + Hqo-s- 

, The first term describes the normal reservoir and reads Hn — X]fe„ a ^kNcl^^Ckj^a, where c|,^^ creates an electron 

' with spin a at the single-particle energy level = efc„ — /iat, with fijy being the chemical potential at the normal 

. reservoir. The coupling to the normal lead is given by the term 

a . 

tN ^ 

which defines F^v = tt Y,kt, l^fe™ " 

All the relevant quantities for the experiment can be extracted from the QD Green's functions in Nambu space 
defined as Gl{t,t') = -ie{t - t'){[l!„[t),^l{t')\^), where = {d„ dl„)^ . In frequency space, the QD Green's 

function can be formally written as G'^{uj)~^ — G^-'^^\iij)~^ — Scr(w), where G^'^°\uj) is the non-interacting dot 
Green's function in Nambu space and the self-energy Eo- takes into account both the coupling to the leads and the 
Coulomb interaction. Of course, the full problem cannot be exactly solved and one needs to resort to approximations. 
In the main text, we present calculations using a Hartree-Fock (HF) decoupling of the self-energy It has been 
demonstrated that such HF decoupling gives reliable results when benchmarked against more sophisticated methods 
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such as numerical renormalization group Q . The HF selfenergy is obtained by considering the first order diagrams in 
the Coulomb interaction. Its diagonal components in Nambu space are given by Sf^^ = — Sf^^^. = [/(n_cr), where 
the spin-resolved occupations {ric) ~ — ^ / <k>jImG\-Y^{bS) have to be calculated self- consistently. Importantly, also 
the anomalous self-energies a — (^21^)* ~ ^^('^t'^i) have to be taken into account. The explicit expression for 
G^^^^(w) reads 

The Andreev level spectrum of the system is given by the roots of the determinant D(uj) = Det\Cr^^^ (lS)^^\., namely 
by the solutions of 

{uj-e^a^ iFw + Fg - t/(n_^))(tj ^ + iV + Fg -t- U {n„)) 

-(rg^^=^ + c/(dtd,))(Fg-^=^ + = 0. (4) 

The QD spectral function is defined as A[uj) = —^ImTr[G]j:^^{uj)], where the trace includes both spin and Nambu 
degrees of freedom. 

In practice, the results presented in the main text were obtained by discretizing the Fourier space in a finite mesh 
of size A'' = 2^^ with oji E [—D, D] and cutoff D = 25A. Starting from an initial effective mean-field solution Sii,-t- = 
—'^22,1 = Y H' iterate the HF equations until good numerical convergence in the spin-resolved occupations is 
reached (as a criterium for good convergence, the iteration stops when the relative error between successive occupations 
in the iteration loop is less than 10^^). 

For a given bias voltage eV = — MSi the conductance across the system is given hy Q = dl /dV, where the total 
current can be decomposed into Andreev and quasiparticle contributions, I — Ia + Iq- For voltages eV < A, the 
quasiparticle current is zero, Iq = 0, and the only contribution comes from Andreev processes. The Andreev current 
reads 

/a = ^ y di^lfNi^ -V)- f^iio + V)]Ta{lo), (5) 
where the Andreev transmission, defined as 

T^H = 4F2^^|Gf2f;Hp, (6) 

describes processes in which an electron (hole) from the normal side is reflected as a hole (electron) while an extra 
Cooper pair is created on the superconducting side. For voltages above the gap, also the quasiparticle contribution is 
finite and is given by: 

Iq^^J MfNioJ -V)- /g(c.)]TQ(w), (7) 
with a quasiparticle transmission defined as 



Tq{uj) = 4r^Fg0(|^| - A) I ^[|Gfi^.(c.)P + IGgU^r - ■^Re{G^,^A^)[G^2U^)n]. (8) 

This quasiparticle contribution consists of three proccesses: 1) conventional tunneling, 2) an electron (hole) in the 
normal side is converted into a hole (electron) excitation in the superconducting side (branch crossing Q), and 3) 
quasiparticle transfer from the normal lead into the superconducting lead, while creating (or annihilating) a Cooper 
pair as an intermediate state. 



3 



Analytic model for energy level repulsion 



Below we derive a simple expression which describes the level repulsion between the doublet states and the states 
in the continuum of the superconducting lead at small coupling <C A, in the regions of gate voltage corresponding 
to a singlet ground state. 

The Hamiltonian for the superconducting lead given by Eq. ([2]) can be diagonalized after a Bogoliubov transforma- 
tion, 

Hs =^J2'^''s"fis<r^ksa, (9) 

where jksa = u^sCksa + (''Vkscl^ks-a ^^'^ annihilation operators of Bogoliubov quasiparticles with energy e^g = 
[A^ + Cfe^l^^^i and Uks, — [(1 ± ^fcg/efcs)/2]^/^ are the BCS coherence factors. The projection of the Hamiltonian 
Hqd~s to a subspace of states with energy close to A in the leads then reads 



Hqd-s = E ( ^ + ^ ) ^L^^^s^ + + 7f E {^L - ^7fcs-.) + h-c] . (10) 



ks 



The sums in the r.h.s of this equation are restricted to momenta such that |^a;sI ^ A, and Uks^Vks ~ l/v^. 

At vanishing tunnel coupling, the ground state at ef > is the product state of the BCS ground state in the lead 
and the empty level in the dot, which we denote |0). Lowest energy excited states consist of the singly occupied level, 
dj.|0), as well as the BCS ground state filled with one Bogoliubov quasiparticle, T^g^l^)- When their energies and 
eks are close and the tunnel couplings Vk^ are finite, the discrete state on the dot and the continuum of states in 
the leads hybridize. Then, we may use 1^'^) = (Adl + J2ks ^fcsTfcso-)!'^) as a variational wavefunction for the excited 
states with spin a. From the set of equations 



{E-e^)A = l=Y.V,*^Bks, (11a) 

(lib) 




that determine possible eigenenergies E associated with the wavefunction l^'a), we obtain the following equation for 
the bound state excitation energy E of state \a): 



- = -Fsv/A/[2(A-C.)]. (12) 



It yields 



'e.~F5v/A/[2(A -£,)], A-e,»(F|A)i/3, 

A-(F|A/2)i/3, ,^^A, (13) 

^ A - F|A/[2(e, - An - A » (F|A)i/3. 

Equation ([T3)) describes the anticrossing (or level repulsion) between the dot state and the BCS continuum. Namely, 
a bound state forms at all values of it gets closer to the edge of the BCS continuum - and the splitting ICi ~ Ctl 
vanishes - as is increased. Note that the excitation energy of the bound state coincides with the Andreev level 
energy obtained from Eq. Q at F5 <C A, in the regions near the edge of the continuum spectrum in the lead, where 
(rif) — (n^) = and the last term in the l.h.s gives a negligible correction in (Fs/A)^/'^ ^ 1. 



Similarly, when e-o- is close to —{U + A), and the dot is doubly occupied in the ground state (i|(i||0) at vanishing 

sAi) 



tunnel coupling, we may use l^'o-) — {Adl- +X]fc ^ksll o-'^t^l)!'^) ^ variational wavefunction at finite coupling and 



obtain the equation 



+ + C/ = -F5v/A/[2(A-C.)]. (14) 
for C,„ E — (2e + U). It is also in correspondence with Eq. ^ in the regions where (ri-|-) — (n^) — 1. 
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Magnetic field dependence of A 



The magnetic field dependence of tlie superconducting gap, A, was estimated from tlie dl/dV vs {B, V) measure- 
ment shown in Suppl. Fig. la, which was taken at the center of an even valley. As discussed in the main text, the 
Andreev levels appear at eV ± A when the system lies deep inside the singlet GS. Thus, A(i3) is readily obtained 
from the evolution of the sub-gap resonances as a function of B (Suppl. Fig. lb). This experimental dependence was 
used for obtaining the theoretical dl/dV plots at finite B (shown in Figs. 2 and 3 of the main text). 




Bx(T) 



FIG. 1: (a) dl /dV vs. {B, V) measurement taken at the center of an even valley, (b) Superconducting gap extracted from (a) 
as a function of Bx- 



Characterization of second S-QD-A'^ device 



Here we present additional data obtained from the device discussed in Figs. 3 and 4 of the main text. The device 
was first characterized in the normal state (Suppl. Fig. 2a). To this end, an external magnetic field B = 2.5 T was 
applied above the critical field i?c ~ 2 T, driving the vanadium electrical leads across the superconducting transition. 
The resulting charge stability diagram revealed a split Kondo resonance within the odd-occupied diamond (highlighted 
by the white arrows). Furthermore, a charging energy U ~ 4.5 meV is extracted from the height of the Coulomb 
diamond, whereas an average tunnel coupling < F > 0.8 meV is estimated from the FWHM of the Coulomb 
resonances. 

The series of dl/dV vs. {Vpg, V) plots shown in Suppl. Fig. 2b [B = 0, 0.5 and 0.75 T) depicts the B-field evolution 
of the sub-gap Andreev levels in the second device. It is noteworthy that the data herein discussed shows a shift of « 
20 meV in the Vpg axis when compared to the corresponding plot in the main text (Fig. 3a). This can be attributed to 
a local charge rearrangement in the vicinity of the QD which, due to capacitive coupling, effectively results in a small 
shift of the charge stability diagram. In spite of the shift, the dl/dV features clearly remain unaltered. The data 
shown in Suppl. Fig. 2 is qualitatively analogous to the situation discussed in the top (experimental) row of Fig. 2 in 
the main text. Specifically, at _B = (left-most panel), the sub-gap resonances cross twice the Fermi level, delimiting 
the boundaries between the singlet {S) and doublet (D) ground states (GSs) . In agreement with the discussion in the 
manuscript, the Andreev levels split with increasing magnetic field when the system lies in the singlet GS. 



Supplementary data on energy level repulsion effect 



Suppl. Fig. 3 contains further information concerning the energy level repulsion effect discussed in the manuscript. 
The right panel of Suppl. Fig. 3a displays an additional dl/dV vs {B, V) measurement taken at the position of the 
dashed line in the left panel, which is further away from the crossing point than position 2 in Fig. 3a of the main 
text. As discussed in the manuscript, the i3-dependence of the Cj, peaks deviates from the expected Zeeman splitting 
of the doublet state, due to the level repulsion of the | l> state by the continuum of quasiparticle states above the 
superconducting gap. Suppl. Fig. 3b demonstrates that this effect becomes more pronounced as the energy of the 
Ci peaks approaches A, as expected from theory. The plotted ^-factor values were estimated from the slopes of the 
vs. B data (red circles and blue triangles, respectively). The horizontal dashed line is positioned at g « 5.5, 
corresponding to the value estimated from the inelastic cotunneling dl/dV steps in the normal state (shown in Suppl. 
Fig. 3c). By its turn, the vertical dashed line signals the position of the singlet-doublet phase boundary. The plot 
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Vpg(mV) Vpg(mV) Vpg(mV) Vpg(mV) 

FIG. 2: (a) Characterization of second device in the normal state {B = 2.5 T). The green hnes are guides to the eye depicting 
the hmits of the Coulomb diamonds with even (E) and odd (O) occupation. The white arrows highlight inelastic cotunneling 
steps in dl /dV related to the split Kondo resonance, (b) Effect of B on the Andreev levels observed in the same device. A 
series of dl /dV vs {Vpg, V) plots taken at _B = 0, 0.5 and 0.75 T is shown. S and D indicate the regions in which the ground 
state is a singlet or a doublet, respectively. 

clearly shows how the g-factor extracted from the Cj, peaks is strongly suppressed when moving away from the crossing 
point. The values obtained from the C-f- peaks, on the other hand, show no significant plunger gate dependence. 




Vpg(mV) Vsd(mV) 

FIG. 3: (a) Right panel: dl/dV vs. [B,V) plot taken at the position marked by the dashed line in the left panel. The 
measurement shown in the left panel is the same as that shown in Fig. 3a in the main text, (b) factor estimated from ^-t- 
and C,i peaks as a function of Vpg. The level repulsion effect is evidenced by the suppressed C,i <?- factor at positions further 
away from the crossing point, where the C,^ peak energy approaches A. (c) Inelastic cotunneling dl /dV steps resulting from the 
field-induced splitting of a spin-1/2 Kondo resonance. This measurement, taken at B = 2.3 T, yields a normal state g ~ 5.5. 



Supplementary data on angular dependence of the Andreev levels at finite B 

To complement the angular dependence data shown in Fig. 4b of the main text, we include here the complete 
_B-field dependences measured at = 0, 7r/2 (Suppl. Fig. 4a). The plot shown in the left panel [9 = 0) is analogous 
to the measurement shown in Fig. 4a of the main text. From this, a QPT field Bgpj, « 0.6 T is estimated. The 
corresponding QPT measurement at ^? = 7r/2 (right panel) shows the same qualitative features, however with two 
significant differences. The most relevant difference is that the slope of the Andreev levels is reduced, indicating a 
smaller g-factor. This results in a QPT field Bgpj, > 1 T. The ij-factor anisotropy is highlighted in Suppl. Fig. 4b, 
where the data points were extracted from the angular dependence of Ct and at a field magnitude \B\ = 0.6 T 
(data shown in Fig. 4b of the main text). The second difference between the two data plots in Suppl. Fig. 4a is in 
the value of the respective critical fields: sa 2T against B^ w 1.5T. 
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(a) dl/dV p^lh] dl/dV [le'lh] 




FIG. 4: (a) _B-field dependences of the Andreev levels measured ai 9 — (left panel) and = 11/2 (right panel). The inset 
in the left panel shows the orientation of the field with respect to the nanowire (6 = corresponds to B aligned parallel to 
NW). (b) g-factor anisotropy measured from the Andreev level splitting at \B\ — 0.6 T. (c) Schematic diagram summarizing 
the angular dependence behaviour. Due to the (;-factor anisotropy, the QPT occurs at different fields for 6 — (solid blue line) 
and 7r/2 (dashed violet line). Lines I and II highlight the dl/dV features observed when \B\ is fixed at |-BqptI s-nd I^qptI) 
respectively, (d) Angular dependence measurement taken with \B\ = \Bqpj, = 0.35 T. Here, the zero-bias QPT peak appears 
for B perpendicular to the nanowire axis. The circles in the line profile plot (right panel) ascribes the measured dl/dV peaks 
to the Andreev levels shown in (c). 



The scheme depicted in Suppl. Fig. 4c summarizes the evolution of the Andreev levels for the longitudinal (solid 
line) and perpendicular (dashed line) field directions, as deduced from the data in panel (a). For a field amplitude 
\B\ — {Bqprpl — 0.6 T, a zero-bias peak (ZBP) is observed at = (or 9 ^ n) in correspondence of the blue dot along 
line I. At 9 = Tr/2 (or 9 = 3tt/2), the ZBP is split into two peaks in correspondence of the violet dots along line I. 
This is indeed the behaviour observed in Fig. 4b of the main text. The scenario is reversed when the field magnitude 
is fixed at \Bqpj,\ (as indicated by line II). Unfortunately, the QPT in the perpendicular direction occurs near 
the closing of the gap (Suppl. Fig. 4a). To better visualize the angular dependence of the Andreev levels in the case 
\B\ = Bqpj,, the Vpg position was moved closer to the singlet-doublet phase boundary, so that the Andreev levels are 
closer to the Fermi level at i? = 0. Suppl. Fig. 4d shows the angular dependence measurement taken at \B\ = 0.35 
T. In this case a ZBP is observed for 9 = Tr/2 (or 9 = 3tt/2), in correspondence of the violet dot along line II, while a 
split peak is observed for = (or 9 — n), as denoted by the pair of blue dots along line II. 
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